Hadron correlators with improved fermions 
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We investigate point-to-point correlation functions for various mesonic and baryonic channels using the O(a)- 
improved Wilson action due to Sheikholeslami and Wohlert. We consider propagators to both time slices and 1. 
We find that discretisation effects are more pronounced than those reported with unimproved Wilson fermions, 
but that the same procedure for removing finite size effects is successful. Extrapolating to the chiral limit, we 
see the notable features predicted phenomenologically: the ratio of interacting to free correlators in the vector 
channel is roughly constant to about 1 fm, while in the pseudoscalar channel the ratio increases markedly due to 
the strong binding. 
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1. Introduction 

Traditionally, much of lattice gauge theory has 
between concerned with calculations of on-shell 
quantities such as masses and matrix elements. 
However, it has been emphasised by Shuryak pi 
that there is a lot of interesting physics to be 
had by looking at off-shell quantities. In the re- 
view cited, the fundamental objects of interest are 
space-like correlations between quantum fields in 
the vacuum: 



K{n-r 2 ) = (0|TJ(ri)J(r 2 )|0), 



1 



where J(r) are currents corresponding to some 
suitable quantum numbers: in meson channels, 
J(r) — ip(r)Tip(r), while in baryon channels, 
J{r) = eijk^p{r)Tiip(r)T2ip(r), for some appropri- 
ate Dirac structure T, T\ and T2. 

These correlations decay very fast; at large vol- 
ume the dominant behaviour is exponential in the 
lightest bound-state mass. It is therefore conve- 
nient to normalise them against the same corre- 
lations calculated for non-interacting quarks, in 
other words using a frozen gauge configuration. 
We therefore define 



R(r) 



K, 



interacting 



(r) 



(2) 



K fTec (r) 

The correlations defined in this fashion contain 
in principle a great deal of information about all 
that is happening in the chosen channel. In this 



way we can show up, for example, the striking 
differences between the vector and pseudoscalar 
which are far from obvious if one looks only at 
masses or decay constants, and hence shed more 
light on the mechanisms underlying confinement 
and chiral symmetry breaking. 

A detailed analysis of point-to-point correla- 
tions has been carried out by Chu et al. 0; we 
refer the reader to that paper for more detailed 
definitions, formulae and methods. Much of our 
methodology is inherited from theirs. Here we de- 
scribe calculations using Sheikholeslami- Wohlert 
0(a)-improved (S-W) fermions Q. In addition, 
our lattice is slightly finer than that of ref. (g), 
though our lattice size is correspondingly smaller. 

We shall now describe our data and its analysis 
and comparisons with phenomenological results, 
then proceed to conclusions. 



2. Analysis 

Our data sets are on 24 3 • 48 quenched lattices 
at inverse coupling /3 = 6.2. The quark prop- 
agators, inverted using the S-W action, are at 
three hopping parameters k = 0.14144, 0.14226 
and 0.14262. On this lattice, K C riticiai = 0.14314. 
We have currently 27 configurations at all three 
k; for consistency we restrict ourselves here to 
these 27. The S-W action adds a two-link part 
to the Wilson action and also requires a rotation 
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Figure 1. Anisotropy effects in the vector channel 



of the fermion fields. This is known to improve 
the perturbative scaling behaviour J4|]; however, 
there is no particular reason to suppose the type 
of analysis we are performing will benefit from 
this. 

We then sum over all rotationally equivalent 
points on all configurations (taking account of an- 
tisymmetry in our baryon correlations), using a 
bootstrap over all data to generate the errors. We 
take the ratio of these with free propagators gen- 
erated using an adaption of the method of Car- 
penter and Baillie || . Initially we used only prop- 
agators from the origin to all sites at t = 0, i.e. 
on the same time slice as the origin. 

We use the method of ref. [pf to improve finite 
size effects by removing images of the source at 
the origin. This works well. 

However, we have no equivalent quantitative 
way of dealing with anisotropy. Examination of 
the data shows that there are strong effects up to 
approximately six lattice units; these appear to be 
significantly larger than those seen by Chu et al. 
with the Wilson action |2J. These authors point 
out that propagators to points nearest the body- 
diagonals of the lattice include more paths and 
should therefore be nearer the continuum values. 
In our case it appears that even these points are 
suffering from significant anisotropy for distances 
r < 6. 



To shed further light on this, we calculate cor- 
relators to points on the next time slice, t = 1. 
We have verified directly that the effect of the dif- 
fering lattice size in the time direction has negligi- 
ble impact on the correlations. We find that these 
points lie consistently lower than the t = points, 
so much so that it is highly unlikely that the 'more 
paths means lower anisotropy' ansatz can be sim- 
ply applied in our case. On the other hand, the 
interpretation 'more paths means lower correla- 
tions' does seem to hold: in figure [l] we show R(r) 
from [n with the image correction applied. The 
different symbols denote co-ordinates with differ- 
ent numbers of zeroes: i.e., (x,y,z,t), (x,y,z,Q), 
etc., with x, y, z and t non-zero. The expecta- 
tion is that the ratio remains roughly constant 
for small r (note that this is for quarks around 
the strange mass, for which the behaviour of the 
anisotropy is cleanest), in other words that from 
r = to 8 the ratio should be ~ f.5. There is 
no easy way to make the data fit this. Presum- 
ably this is an effect primarily of the rotation of 
the fermion fields. Indeed, the major additional 
problem appears to be that our correlations plum- 
met for some small values of r; apart from that 
the anisotropy is similar to the Wilson case. 

Matters are improved slightly by the choice of 
Chu et al. to restrict points to those making an 
angle 8 with the body diagonal of the spatial lat- 
tice with cos(6>) > 0.9. In this case the scatter 
of points is consistent with the statistical scatter 
from about r ~ 6, and our fits are performed from 
this distance on. 

2.1. Chiral extrapolations 

We follow the usual UKQCD convention in ex- 
trapolating our results to zero quark mass, rather 
than the alternative of extrapolating to the phys- 
ical pion mass. Our chiral extrapolations are lim- 
ited by the fact that currently we only have three 
quark masses. We have tried various forms of chi- 
ral extrapolation and have settled on a fit which 
linearly extrapolates log(i?(r)). 

2.2. Fits to phenomenological forms 

A reasonable phenomenological approximation 
to the data should be found be assuming a smooth 
continuum background at large energy plus some 



number of resonance poles. In practice one pole is 
entirely adequate and it is unlikely that our data 
could be persuaded to yield a convincing second 
pole. For the channels shown here this gives four 
parameters to fit: an overall normalisation, the 
position (m) and height (A) of the delta function 
pole, and the threshhold for the continuum. Be- 
cause we cannot fit for r < 0.6 the estimation of 
the continuum threshhold is generally poor and 
the fits shown just use a reasonable fixed value. 
Details of the forms are as in ref. M . Errors are 
found by jackknife, completely re-analysing the 
data with each set missing in turn. 

We show the extrapolated data with the 
cos(9) > 0.9 cut off and fits for the vector, nu- 
cleoli and pseudoscalar in figure @. One fermi 
corresponds to about fourteen of our lattice spac- 
ings. We find a slightly low rho mass (620 ± 90 
MeV) and a high nucleon mass (1170 ± 70 MeV); 
also the pseudoscalar is rather heavy (300 ± 20 
MeV), throwing some doubt on the validity of 
extrapolating to zero quark mass for this type of 
analysis. 
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3. Conclusions 

We have examined point-to-point hadron cor- 
relation functions using an O(a) improved action. 
We find more severe finite lattice spacing effects 
than reported using Wilson fermions. Nonethe- 
less, the physical features of the results are sim- 
ilar: asymptotic freedom in the vector channel 
persists to about 1 fm, but is violated quickly in 
the pseudoscalar channel. 
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Figure 2. Extrapolated data and fits for the pseu- 
doscalar, vector and nucleon. The vertical scale 
is arbitrary: anisotropy makes normalisation dif- 
ficult. 



